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Abstract 

Recent measurements of very low-energy (pl < 100 MeV/c) p annihilation 
on light nuclei reveal apparent suppression of annihilation upon increasing 
the atomic charge Z and mass number A. Using p-nucleus optical potentials 
V pt, fitted to p-atom energy-shifts and -widths, we resolve this suppression 
as due to the strong effective repulsion produced by the very absorptive V op t ■ 
The low-energy p-nucleus wavefunction is kept substantially outside the nu- 
clear surface and the resulting reaction cross section saturates as function of 
the strength of Im V op t- This feature, for E > 0, parallels the recent predic- 
tion, for E < 0, that the level widths of p atoms saturate and, hence, that 
p deeply bound atomic states are relatively narrow. Antiproton annihilation 
cross sections are calculated at pi = 57 MeV/c across the periodic table, and 
their dependence on Z and A is classified and discussed with respect to the 
Coulomb focussing effect at very low energies. 
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I. INTRODUCTION 



Experimental results for antiproton annihilation cross sections at very low energies (pi < 
100 MeV/c), below the pp —>■ fin charge-exchange threshold, have recently been reported 
for light nuclei [0-0]. At these energies the total p reaction cross section consists only of p 
annihilation. Bianconi et al. || reported the first ever measurements of p annihilation on 
Ne in the momentum range of 53-63 MeV/c. Comparing with other data, they showed that 
whereas at relatively higher energies (p^ ~ 200 — 600 MeV/c) ratios of p annihilation cross 
sections on different nuclei exhibit the well-known A 2 ! 3 strong-absorption dependence, such 
ratios at very low energies defy any simple, obvious regularity. It has been demonstrated 
that the 'expected' ZA 1 ^ 3 dependence of these cross sections on the atomic charge Z and 
mass number A is badly violated ]3H. Antiproton annihilation cross sections at very low 
energies simply do not rise with A as fast as is anticipated. For example, over the whole 
mass range studied so far (H, D, 4 He, Ne), the Ne/H ratio of total annihilation cross sections 
is at least 6 times smaller than expected. 

In the present Letter we study low-energy p annihilation on nuclei, using the opti- 
cal model approach. Optical potentials have been very successful in describing strong- 
interaction effects in hadronic atoms || , including p atoms || . It has been noted, for pions, 
that the total reaction cross sections at low energies are directly related to the atomic-state 
widths, and that once a suitable optical potential is constructed by reasonably fitting it to 
the atomic level shifts and widths in the negative energy bound-state domain, these total 
reaction cross sections are reliably calculable [|]||. The recent publications |]-[3|] of experi- 
mental results of total cross sections for p annihilation on nuclei at very low energies raise 
the intriguing possibility of connecting these two energy regimes in a systematic way also for 
antiprotons. However, most of the data on annihilation cross sections for p are for very light 
nuclei, where the concept of a rather universal optical potential that depends on A and Z 
only through the nuclear densities is questionable. For this reason we use optical potentials, 
in the present work, mostly for crossing the E = borderline within the same atomic mass 
(A) range, from bound states to scattering. This optical model approach is different from 
the allegedly model-independent scattering length approximation || which we have found 
to be less useful. We defer a discussion of the latter statement to a separate publication. 

Among the targets used in the p annihilation measurements PJ-BJ, the only 'real' nuclei 
are 4 He and Ne, in the sense that the conventional optical model is not expected to be applied 
to the lighter targets of hydrogen or deuterium. In the present work we show that the very 
low-energy p total annihilation cross sections on 4 He and on Ne || are reproduced well 
by reasonable optical potentials fitted separately to the p atomic data at the relevant atomic 
mass range. These potentials are strongly absorptive, which leads to a remarkable saturation 
of the total reaction cross section with increasing A. Strong absorption has very recently 
been shown |TU| , |TT| to lead also to saturation of the widths of p atomic states and to the 
prediction of relatively narrow deeply bound p atomic states. In view of the close analogy 
between bound-state widths and total reaction cross sections (see Eqs. (JTJ.0) below), it should 
come as no surprise that the underlying physics is the same. In fact, several pD calculations 
12] {L4J have found that the p absorptivity in deuterium, as suggested by the magnitude 



of the imaginary part of the pD s-wave scattering length, is weaker than in hydrogen due 
partly to the repulsive effect of absorption. Lastly, we use p optical potentials fitted to a 
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comprehensive set of atomic data across the periodic table to predict p total annihilation 
cross sections at very low energies for a wide range of A values, with the exception of very 
light targets. We show such calculated cross sections at pi = 57 MeV/c and discuss their A 
and Z dependence, particularly with respect to the ZA 1 ^ dependence arising from Coulomb 
focussing in the limit of very low energies. We also demonstrate the extent to which this 
dependence is violated at 57 MeV/c. 



II. OPTICAL POTENTIALS 

In the present work we aim at connecting p atoms, at energies slightly below threshold, 
with p annihilation on nuclei at very low energies above threshold, using an optical potential 
V opt . Assuming for simplicity a Schrodinger-type equation, the width T of an atomic level 
is given (non-perturbatively) by: 

r _ JImK P t(r)|^(r)| 2 rfr 

2 " /|^(r)| 2 dr ' 1 ' 

Here ^(r) is the p full atomic wavefunction. The corresponding expression for the total 
reaction cross section at positive energies is 

<tr = ~J |x(r)| 2 ImK P t(r)rfr . (2) 

Here x( r ) is the p - nucleus elastic scattering wavefunction and v is the cm. velocity. 

There exist several optical model potentials |||| parameterized in terms of nuclear den- 
sities which are quite successful in reproducing the strong interaction effects in not too light 
nuclei, e.g. fitting all data for nuclei heavier than carbon. However, most of the good quality 
data on p annihilation at very low energies are for very light nuclei where a straightforward 
application of the potential deduced from heavier targets is not expected to yield good fits. 
We have therefore started the present analysis by studying p atoms of 3 ' 4 He. 

The p nucleus optical potential used here is given by the Hp' expression [H 

2fj,V opt (r) = -4tt(1 + £)Mr) , (3) 

where m is the nucleon mass, \l is the reduced p-nucleus mass, b is a complex parameter 
obtained from fits to the data and p(r) is the nuclear density distribution normalized to A. 
Trying to fit level shift and width data in 3 ' 4 He [15| for the 2p and 3c? states simultaneously, 
we always ended up with the calculated width of the 2p state in 4 He being too small, thus 
contributing an unacceptably large value to the total x 2 of the fits. Handling each isotope 
separately did not change this situation and we then decided to fit only the data for the 
2p states for the two isotopes. This is somewhat unfortunate as one expects to be able to 
use /-independent potentials, as is indeed the case with heavier targets [fHI . However, since 
our main concern here was to use atomic potentials at (very low) positive energies where d 
and higher partial waves contribute very little in light nuclei, we considered this procedure 
acceptable. Note that we are using the same potential for all the partial waves in the positive 
energy regime. 
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Very good fits to the data with reasonable values for the complex parameter 60 could 
be obtained only when a 'finite range' was introduced in the form of a Gaussian folding of 
a pN interaction into the nuclear density distributions. This is in contrast to the case of 
heavier targets where such a procedure was not necessary. Nevertheless, we consider the 
resulting potential, with 6o=0.49+z3.0 fm and a Gaussian with a range parameter of 1.4 fm 
folded into the nuclear density distribution (hereafter referred to as potential (a)), as quite 
reasonable. It fits within the errors the measured strong interaction shifts and widths of the 
2p state in antiprotonic 3 He and 4 He. 



III. RESULTS AND DISCUSSION 

Turning to positive energies, we first calculate the total reaction (annihilation) cross 
section for 57 MeV/c p on 4 He, in order to make it possible to compare results for He and 
Ne [Q. Using the above potential (a), the calculated cross section is 901 mb, in excellent 
agreement with the interpolated value of 915±39 mb 0. For incident momenta of 47.0 and 
70.4 MeV/c we calculate values of 1060 and 771 mb, respectively, where the experimental 
values are 979±145 and 827±38 mb, respectively. Reasonably good agreement is therefore 
established between experiment and predictions made with a potential derived from fits to 
p atoms of 3 ' 4 He. 

Figure [l] demonstrates the extreme strong-absorption conditions which are relevant to 
the p nucleus interaction at very low energies (and for p atoms). It shows calculated reaction 
cross sections for p at 57 MeV/c on 4 He and Ne as function of the strength Im bo of the 
imaginary part of the potential (a) described above, with the rightmost edge corresponding 
to its nominal value Im bo=3.0 fm. It is seen that as long as the absorptivity (Im bo) is 
very weak, less than 1% of its nominal value, <jr is approximately linear in Im bo, which 
according to Eq. (^) means that the p wavefunction depends weakly on Im bo- However, 
already at below 5% of the nominal value of Im bo the reaction cross sections begin to 



saturate, much the same as for the widths of deeply bound p atomic states [10,11|. The 
mechanism is the same in both cases, namely exclusion of the wavefunction from the nucleus 
due to the absorption, which reduces dramatically the overlap with the imaginary potential 
and consequently reduces the integrals in Eqs. (0,0). The onset of saturation is determined 
approximately by the strength parameter 2/z(Im V^ pt )-R 2 , where R is the radius of the nucleus. 
Thus, saturation of ctr in Ne starts at a smaller absorptivity than its onset in 4 He. The 
Ne/ 4 He ratio of ctr values changes, due to this effect, from about 15 in the perturbative 
regime of very weak absorptivity to about 3 in the strong-absorption regime. However, the 
precise, detailed pattern of the change also depends on the real part of the optical potential 
which may increase or reduce the exclusion of the p wavefunction from the nucleus. 

Figure |] shows calculated p reaction cross sections at 57 MeV/c across the periodic table. 
The dot-dashed line is for the above mentioned potential (a) which is expected to be valid 
only in the immediate vicinity of He. The dashed line (b) is for the first potential from Table 
7 of Ref. H which fits p atom data over the whole periodic table, starting with carbon. This 
potential is not expected to fit data for very light nuclei, and indeed it does not fit the 4 He 
annihilation cross section. However, it is noteworthy that for A > 20 the two potentials 
predict almost the same cross sections and certainly display a very smooth and similar 
dependence on A. A smooth A dependence of the calculated cross sections is generally 
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expected once several partial waves contribute, as we elaborate below. Also shown in the 
figure is the recent experimental result || for Ne, with very limited accuracy. Furthermore, 
for A > 20 the points along the solid line are the calculated n - nucleus total reaction cross 
sections, obtained from potential (b) by switching off the Coulomb interaction. The solid 
line is a fit to an A 1 ^ 3 power law which appears to be appropriate to strong absorption of 
uncharged particles at these very low energies, in contrast to the prediction of approximate 
A independence made in Ref. [[nj for n annihilation at ultra low energies where only s waves 
contribute. We checked that by reducing the absorptivity by 3 - 4 orders of magnitude, a 
proportionality of ctr to A as expected for weak absorption, is indeed observed. Comparing 
the dashed line for negatively charged particles with the solid line for uncharged particles, it 
is clear that the <jr values obtained by including the Coulomb interaction are considerably 
enhanced with respect to the ones obtained without it. Attempting to fit the enhancement by 
a Z a power law, we find that a varies roughly between 1/3 for Ne to 1/2 for Pb. Attempts 
to identify more precisely this enhancement across the periodic table have been recently 
addressed in Ref. fL7| 



Finally, we discuss the ZA 1 / 3 scaling considered in Refs. pfj] and in references quoted 
therein. We note that an attractive Coulomb potential causes focussing of partial-wave 
trajectories onto the nucleus, an effect which at very low energies may be evaluated semi- 



classically (see Ref. |18| for the applicability of the semiclassical approximation to Coulomb 
scattering). The maximal orbital angular momentum Z max , for which the p- nucleus distance 
of closest approach is smaller than the nuclear radius R, at sufficiently low energy is given 
by 

(Z max + 1/2) 2 « 2 V kR , (4) 
where k is the cm. momentum and 

r) = 2Z/k L a ( i ) (5) 

is the corresponding Coulomb parameter. Here k^ = pt/%, where pi is the p laboratory 
momentum, and is the Bohr radius of the pp atom. Since 2r] >> kR at very low energies 
and for high values of Z, the value of Z max greatly exceeds the value kR which is appropriate 
to uncharged particles (77 = 0). For example, on Pb at pl = 57 MeV/c, Z max is about 6 for p, 
but is less than 2 for n. Therefore, for strong absorption, where all the lower partial waves 
are totally absorbed, one obtains 

7T ' max 2/7 

ff R~pE( 2; + 1 )~^» rf > (6) 

which we rewrite as 

a R n (Anro/a^iZA^/kkL) , (7) 

with ro defined by R = tqA 1 ^ 3 . A similar approach was used by Faldt and Pilkuhn [ 19 



to calculate Coulomb corrections to ^-nucleus total cross sections near the 3,3 resonance. 
In the present case of very low energies their correction becomes the dominant term, and 
the resulting reaction cross section on the r.h.s. of Eq. @ is substantially higher than the 
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black-disk value <jr = ttR 2 which for uncharged projectiles may be perceived as the unitarity 
limit of gr. Indeed, the calculated saturation plateau values of gr in Fig. [l] are several times 
larger than the appropriate black-disk values. 

In Figure § we plot, for potential (b), ratios of calculated total reaction cross sections at 
57 MeV/c to the scaling parameter S = Z A 1 ! 3 Ikk^. It is clear that gr does not scale with 
S across the periodic table. Also plotted is the asymptotic value 

a R /S -+ Airro/a® « 0.273 , (8) 

with tq = 1.25 fm and = 57.6 fm. For a given p beam energy this asymptotic value, 
arising from Coulomb focussing, is approached very slowly upon increasing Z. 

Before concluding we briefly mention other derivations of the ZA 1 ! 3 dependence of gr. 
The proportionality to Z is usually derived from the well-known Coulomb enhancement 
Gamow factor which at low energies is well approximated by 2-kt] [|l7j]. However, at these 



very low energies, the Coulomb wavefunctions are far from being constant over the nuclear 
volume, so that the applicability of the Gamow factor appears questionable. The origin of a 
possible proportionality of gr to A 1 / 3 is even less clear: for s-wave dominated p annihilation, 
as considered in Refs. |||]] , it is due to the imaginary part of the p scattering length which is 
assumed to closely follow the nuclear radius R. However, for p strongly absorptive potentials 
this imaginary part has invariably been found to be about 1 fm across the periodic table 
PD] (see also Ref. |]]), considerably less than R. It is inconceivable that a proportionality 
to R can arise unless several partial waves contribute significantly. 



IV. CONCLUSIONS 

In conclusion, we have shown that the recently reported annihilation cross sections for 
p on 4 He and Ne at 57 MeV/c are reproduced very well by optical potentials which fit p 
atomic data in the respective mass ranges (Fig. ^j. For these strongly absorptive potentials, 
the asymptotic ZA 1 / 3 dependence of the total reaction cross section was derived from the 
Coulomb focussing effect at very low energies, but at 57 MeV/c it was found not to be 
satisfied across the periodic table except for unrealistically high values of A (Fig. |3|). The 
correct dependence is determined by the saturation property of total reaction cross sections 
at very low energies (Fig. [I]). The apparent suppression of p total annihilation cross sections 
measured at p^ < 100 MeV/c is a direct consequence of the strong absorption, and in 
particular its saturation property. Predictions have been given for p total annihilation cross 
sections at 57 MeV/c across the periodic table which would complement and enrich the 
information deduced from the measurements already available for p atoms. 

A more extended discussion of optical model calculations and their implications for p 
atoms and p annihilation on nuclei, including a connection to pp annihilation at p^ < 200 
MeV/c, is included in a forthcoming publication [ETI . 



CJB wishes to thank the Hebrew University for support for a visit during which this 
work was started. 

This research was partially supported by the Israel Science Foundation. 



5 



REFERENCES 



A. Zenoni et al, Phys. Lett. B 461 (1999) 405. 
A. Zenoni et al., Phys. Lett. B 461 (1999) 413. 
A. Bianconi et al, Phys. Lett. B 481 (2000) 194. 

A. Bianconi, G. Bonomi, M.P. Bussa, E. Lodi Rizzini, L. Venturelli, A. Zenoni, Phys. 
Lett. B 483 (2000) 353. 

C.J. Batty, E. Friedman, A. Gal, Phys. Rep. 287 (1997) 385. 
C.J. Batty, E. Friedman, A. Gal, Nucl. Phys. A 592 (1995) 487. 
K. Strieker, H. McManus, J.A. Carr, Phys. Rev. C 19 (1979) 929. 
E. Friedman et al., Phys. Lett. B 257 (1991) 17. 

K.V. Protasov, G. Bonomi, E. Lodi Rizzini, A. Zenoni, Eur. Phys. J. A 7 (2000) 429. 

E. Friedman, A. Gal, Phys. Lett. B 459 (1999) 43. 

E. Friedman, A. Gal, Nucl. Phys. A 658 (1999) 345. 

S. Wycech, A.M. Green, J.A. Niskanen, Phys. Lett. 152 B (1985) 308. 

CP Latta, PC. Tandy, Phys. Rev. C 42 (1990) R1207. 

V.A. Karmanov, K.V. Protasov, A. Yu. Voronin, Eur. Phys. J. A (2000), in press ( nucl- 
th/0006041| ). 



M. Schneider et al., Z. Phys. A 338 (1991) 217. 

Ye.S. Golubeva, L.A. Kondratyuk, Nucl. Phys. B (Proc. Suppl.) 56A (1997) 103. 

A. Bianconi, G. Bonomi, E. Lodi Rizzini, L. Venturelli, A. Zenoni, Phys. Rev. C 62 

(2000) 014611. 

L.I. Schiff, Quantum Mechanics, third edition (McGraw Hill, New York, 1968) 
G. Faldt, H. Pilkuhn, Phys. Lett. 40 B (1972) 613. 
C.J. Batty, Nucl. Phys. A 411 (1983) 399. 

C.J. Batty, E. Friedman, A. Gal, to be submitted to Nucl. Phys. A (2000). 



6 



10000 



FIGURES 




r| I i i i i i i i i 

0.001 0.010 0.100 1.000 10.000 

Im b (fm) 

FIG. 1. Calculated total reaction cross sections for 57 MeV/c p on 4 He and Ne as function of 
the strength parameter Im 6q of the optical potential (a). 
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FIG. 2. Calculated p total reaction cross sections (open circles) at 57 MeV/c for potentials (a) 
and (b), and for potential (b) but without the Coulomb interaction. Also shown are the two data 
points for 4 He and Ne. 
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FIG. 3. Ratios of calculated p total reaction cross sections at 57 MeV/c, for potential (b), to 
the scale function S. The dashed straight line provides the asymptotic value expected for ctr/S, 
see text. 
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